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Abstract 

The  theory  of  thermoelastic  materials  undergoing  solid-solid  phase 
transformations  requires  constitutive  information  that  governs  the 
evolution  of  a  phase  boundary.  This  is  known  as  a  kinetic  relation 
which  relates  a  driving  traction  to  the  speed  of  propagation  of  a  phase 
boundary.  The  kinetic  relation  is  prescribed  in  the  theory  from  the 
onset.  Here,  though,  a  special  kinetic  relation  is  derived  from  an 
augmented  theory  that  includes  viscous,  strain  gradient  and  heat  con¬ 
duction  effects.  Based  on  a  special  class  of  solutions,  namely  traveling 
waves,  the  kinetic  relation  is  inherited  from  the  augmented  theory  as 
the  viscosity,  strain  gradient  and  heat  conductivity  are  removed  by  a 
suitable  limit  process. 
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1  Introduction 


Materials  can  exist  in  different  solid  phases  which  are  characterized,  at  a 
microscopic  level,  by  their  underlying  crystalline  structure.  A  change  in 
temperature,  stress  or  other  variables  can  induce  a  material  to  change  from 
one  solid  phase  to  another.  Phase  changes  are  usually  classified  according  to 
the  main  mechanism  driving  the  transformation,  viz.,  by  diffusion  of  chemical 
species  or  by  displacement,  although  these  are  not  mutually  exclusive  (see 
Christian  [8]).  The  present  work  deals  with  the  latter  type  in  which  the 
transition  from  one  phase  to  another  is  characterized  by  a  sudden  change  in 
the  crystalline  lattice,  achieved  without  diffusion  but  rather  via  a  cooperative 
movement  of  atoms. 

Inside  a  body,  the  interface  between  two  distinct  phases,  known  as  a 
phase  boundary,  can  be  modeled  as  a  surface  of  zero  thickness  across  which 
the  displacement  is  continuous  but  the  deformation  gradient  suffers  a  jump. 
Other  models,  related  to  the  Landau-Ginzburg  theory,  consider  the  interface 
as  a  transition  layer  of  small  but  finite  thickness  where  the  deformation  gra¬ 
dient  varies  rapidly  but  continuously  (see,  e.g.,  Ericksen  [13],  PENROSE  & 
Fife  [21],  Truskinovsky  [27]).  Implicit  in  the  term  “small”  is  the  idea 
that  the  model  has  some  kind  of  length  scale.  Some  related  theories  that 
incorporate  separate  balance  principles  associated  with  the  motion  of  mate¬ 
rial  particles  and  the  motion  of  the  phase  boundary  have  also  been  proposed 
(see,  Fried  &:  Gurtin  [15],  Lusk  [19]). 

Within  the  framework  of  classical  thermoelasticity  theory,  the  thermo¬ 
mechanical  properties  of  a  material  can  be  described  by  a  Helmholtz  free 
energy  density  that  depends  smoothly  on  the  deformation  gradient  and  the 
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temperature.  In  order  to  model  a  material  that  can  exist  in  different  solid 
phases — and  hence  analyze  phase  transitions — the  Helmholtz  potential,  at 
some  temperatures,  exhibits  multiple  wells.  Each  well  is  associated  to  a  dis¬ 
tinct  solid  phase.  The  nonconvexity  of  the  potential  allows  the  existence  of 
weak  solutions  to  dynamic  problems.  The  solutions  involve  shock  waves — 
across  which  the  deformation  gradient  is  discontinuous  but  the  material  on 
each  side  of  the  surface  of  discontinuity  is  in  the  same  phase — and  propagat¬ 
ing  phase  boundaries  if  on  each  side  the  material  phases  are  different. 

Phase  transformations  in  solids  have  been  studied  in  recent  papers  by 
Abeyaratne  Knowles  [1,  3,  4,  5].  In  particular,  [5]  deals  with  the 
adiabatic  theory  for  a  propagating  phase  boundary  in  a  one-dimensional 
thermoelastic  solid  where  an  initial  value  problem  of  the  Riemann  type  is 
analyzed.  The  phase  boundary  is  modeled  as  a  sharp  interface.  It  is  found 
that  the  solution  to  the  Riemann  problem  is  not  unique.  In  order  to  recover 
uniqueness,  Abeyaratne  and  Knowles  propose  the  existence  of  an  additional 
piece  of  constitutive  information,  a  kinetic  relation,  which  restricts  the  mo¬ 
bility  of  the  phase  boundary  and  singles  out  a  unique  solution  for  the  Rie¬ 
mann  problem.  Following  Onsager’s  formalism  in  thermodynamics  (see,  e.g., 
Callen  [6],  Truesdell  [26]),  the  kinetic  relation  relates  a  driving  traction 
(associated  with  a  jump  in  entropy)  to  the  speed  of  propagation  of  the  phase 
boundary  and  to  a  state  parameter,  namely  the  temperature.  In  the  sharp 
interface  theory,  the  kinetic  relation  is  introduced  from  the  onset. 

The  purpose  of  the  present  analysis  is  to  derive  a  special  kinetic  relation 
from  an  admissibility  criterion  related  to  a  traveling  wave  problem  in  which 
an  augmented  theory  is  used.  In  the  augmented  (or  regularized)  theory,  the 
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sharp  discontinuities  are  replaced  by  transition  layers  of  finite  thickness.  By 
introducing  dissipation  and  dispersion  mechanisms  associated  with  higher 
order  derivatives  of  the  strain  (viscosity  and  strain  gradient  in  this  case),  the 
phase  boundary  acquires  a  structure,  i.e.,  it  is  possible  to  simulate  a  region 
where  the  strain  varies  rapidly  but  continuously.  The  augmented  theory  also 
includes  heat  conduction  in  order  to  have  a  continuously  varying  temperature 
across  the  transition  layer.  In  a  sense,  the  sharp  interface  theory  should 
be  obtained  from  the  augmented  theory  through  a  suitable  limit  process 
where  the  dissipation  and  dispersion  mechanisms  and  the  heat  conduction 
are  removed. 

In  Section  2,  the  basic  field  equations  and  jump  conditions  are  derived  for 
the  sharp  interface  theory  and  the  results  are  specialized  using  a  model  for  a 
thermoelastic  material  proposed  in  [3].  In  Section  3,  an  augmented  theory  is 
introduced.  Some  general  features  of  solutions  to  the  traveling  wave  problem 
are  derived  in  Section  4.  It  is  also  shown  that  the  admissibility  criterion  for 
the  existence  of  a  subsonic  traveling  wave  in  the  augmented  theory  provides 
an  additional  restriction  on  the  data  that  can  be  cast  as  a  kinetic  relation. 
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2  Sharp  interface  theory:  adiabatic  case 

2.1  Beisic  equations. 

Consider  a  bar  that  occupies  the  interval  [— L,  L]  in  a  reference  configuration. 
Consider  longitudinal  motions  of  the  bar  during  a  time  interval  [to,  ti]-  Points 
in  the  reference  configuration  are  denoted  by  {x,  t)  G  [—L,  L]  x  [to,  ti].  Assume 
that  the  particle  x  is  mapped  to  the  position  y{x,  t)  at  time  t,  (i.e.,  y{x,  t)  = 
x+u{x,  t),  where  u  is  the  displacement);  assume  that  u  €  C®([— L,  L]  x  [to,  ti]) 
and  piecewise  C‘^{\-L,L]  x  [to,ti]).  Let  y  =  and  v  —  Ut  he  the  strain 
and  particle  velocity  (the  subscript  refers  to  partial  differentiation).  The 
restriction  -1  <  y{x,t)  ,  Vx,t,  guarantees  that  the  deformation  y{x,t)  is  one- 
to-one.  Let  p  be  the  referential  mass  density  (assumed  to  be  independent 
of  x),  a  the  stress,  e  the  internal  energy  per  unit  mass  and  rj  the  entropy. 
Assume  that  a,  e  and  y  are  piecewise  C^([—L,L]  x  [foi^i])-  The  balance 
of  linear  momentum,  balance  of  energy  and  the  Clausius-Duhem  inequality 
(dissipation  inequality)  for  an  adiabatic  thermomechanical  process  in  a  one¬ 
dimensional  bar  are 


r(f)  =  ^J  0 

Vf  G  [to,ti]  and  V[xi,X2]  C  [-L,  L].  Here,  T{t)  is  the  entropy  production  rate 
in  the  interval  [xi,X2].  The  local  versions  of  (1),  (2)  and  (3)  at  points  where 
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<7,  V,  7,  e  and  rj  axe  smooth  are 


II 

1 

(4) 

{av)^  =  , 

(5) 

r)t>0  . 

(6) 

The  compatibility  equation  is 

Vx  =  It  ■ 

(7) 

Equation  (5),  with  the  use  of  (4),  can  be  expressed  as 

<yit  =  P^t  ■ 

(8) 

Let  X  =  s{t)  be  a  point  in  [xi,  xa]  across  which  some  (or  all)  of  the  fields  are 
discontinuous.  Localization  at  x  =  s{t)  of  the  compatibility  equation  (7)  and 
the  global  balances  (l)-(3)  provide  the  corresponding  jump  conditions,  viz., 


[vj  +  il7l  =  0  , 

(9) 

[o-J  +  psful  =  0  , 

(10) 

[cruj  +  ps|e  +  =  0  , 

(11) 

lv]s  <  0  , 

(12) 

where  s  =  ds/dt  is  the  speed  of  propagation  of  the  discontinuity  and,  for 
any  function  g,  |^]  =  -  g~  and  g^  =  lima;_,.s(t)±  g{x,t).  The  superscript 
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refers  to  one  side  of  the  moving  discontinuity  as  follows:  the  side 

is  the  front  (back)  state  if  s  >  0  or  the  back  (front)  state  if  s  <  0.  Moving 
discontinuities  are  said  to  be  compressive  if  the  strain  in  the  front  state  is 
higher  than  in  the  back  state  and  expansive  otherwise.  The  velocity  jump 
can  be  eliminated  from  (10)  by  using  (9),  i.e., 

[a]  =  .  (13) 

Similarly,  (11),  with  the  use  of  (9)  and  (10),  can  also  be  expressed  as 

(p|el  -  (o->|7])s  =  0  ,  (14) 

where,  for  any  function  g,  (g)  =  (l/2)(p'^  +  g~). 

2.2  Thermoelastic  material. 

A  thermoelastic  material  can  be  characterized  by  its  Helmholtz  free  energy, 
which  is  a  function  of  strain  and  temperature,  or  by  its  internal  energy,  which 
is  a  function  of  strain  and  entropy.  It  can  also  be  characterized  by  its  Gibb’s 
energy  or  by  the  enthalpy.  One  of  these  potentials  is  taken  as  a  fundamental 
quantity  and  the  others  can  be  obtained  from  it  via  various  Legendre  trans¬ 
formations.  If  certain  invertibility  requirements  between  variables  are  met, 
it  is  possible  to  express  these  potentials  in  terms  of  any  combination  of  two 
variables.  If  the  strain  and  the  temperature  are  used  as  variables,  then  the 
stress  and  the  entropy  are  viewed  as  functions  of  7  and  6.  As  shown  in  sub¬ 
sequent  sections,  for  the  thermoelastic  materials  considered  here — which  are 
capable  of  phase  transformations — the  stress  a  is  not  a  monotonic  function 
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of  7  or  0  and  the  entropy  f)  is  not  a  monotonic  function  of  7.  Nevertheless, 
it  is  assumed  that  the  entropy  is  a  monotonic  function  of  the  temperature, 
hence  either  (7, 6)  or  (7, 7/)  can  be  used  as  variables.  The  following  notation 
is  employed:  a  function  of  strain  and  temperature  is  denoted  as  g  and  a  func¬ 
tion  of  strain  and  entropy  is  denoted  as  g.  Throughout  the  present  analysis, 
the  preferred  variables  are  7  and  9.  Let  the  Helmholtz  potential  be  given  by 
'tp  =  9).  This  potential  is  related  to  the  internal  energy  through 


Ip  =  e-9T)  .  (15) 

For  a  classical  thermoelastic  material,  the  stress  a  and  the  entropy  7/  are 
given  by 


0-  =  <7(7,0)  =  pip^i'y,  9)  , 

(16) 

V  =  0(7, 0)  =  -Mk,  0) . 

(17) 

The  isothermal  elastic  modulus  p,  the  specific  heat  at  constant  strain  c  and 
the  coefficient  of  thermal  expansion  a  are  defined  by 

p  =  pip^^,  c  =  -9ip6e,  d  =  -^.  (18) 

V'li 
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Alternatively,  one  can  introduce  the  modified^  Griineisen  coefficient  defined, 
in  terms  of  the  above  quantities,  as 

G  =  ^  ^ 

9ipee  pc 

As  a  fundamental  assumption,  only  materials  with  positive  specific  heat  at 
constant  strain  and  positive  coefficient  of  thermal  expansion  axe  considered, 
i.e.,  V(7, 0)  c{‘y,6)  >  0.  Hence,  it  is  assumed  that  G{'y,9)  >  0.  Moreover, 
since  0  >  0  and  by  (18)2,  then  77  is  a  strictly  increasing  function  of  9.  Thus, 
it  is  possible  to  define  an  inverse  function  9  —  9(^,7]).  The  natural  variables 
for  the  internal  energy  are  7  and  rj.  Prom  equation  (15),  in  terms  of  the 
Helmholtz  potential,  e  is  given  by 

Kit  n)  =  s(t.  ri))  +  e(i,  ri)ji .  (20) 


(19) 


Therefore,  from  (16)  and  (20),  the  stress  can  also  be  expressed  as  a  function 
of  7  and  rj,  i.e.,  a  =  a{'^,ri)  =  ptyipf^r}).  Similarly,  the  temperature  is 
given  hy  9  =  9{'),t])  =  €^(7,7y).  The  isentropic  elastic  modulus  /Xe  is  defined 
as  /2e  =  pc-,-^-  Henceforth,  though,  from  the  invertibility  of  rj  and  9,  the 
internal  energy  e  is  to  be  considered  as  given  by  a  function  e  of  7  and  9,  i.e., 
e  =  e(7, 9)  =  e(7, 77(7, 9)).  In  terms  of  7  and  9,  the  isentropic  elastic  modulus 
He  can  be  expressed  as 


He  P 


i’ee 


(21) 


^The  Griineisen  coefficient  is  usually  defined  as  (1  +  'y)G,  where  G  is  given  by  (19).  It 
turns  out  that  (19)  is  a  more  convenient  parameter.  See  Clifton  [9]. 
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The  isothermal  and  the  isentropic  sound  speeds,  a  and  Oe,  are  defined,  when 
/x,/ie  >  0,  by 

(22) 

Using  the  material  parameters  defined  by  (18),  it  follows  from  (21)  and  (22) 
that  the  sound  speeds  are  related  by 

al  =  a^(l  +  ^ej  ;  (23) 

hence,  if  /i,/ie  >  0  and  since  c  >  0,  then  >  a^.  In  the  adiabatic  the¬ 
ory,  discontinuities  are  classified  into  two  types,  viz.,  shock  waves  and  phase 
boundaries.  Shock  waves  are  related  to  a  discontinuity  where  the  material  on 
each  side  is  in  the  same  phase.  A  phase  boundary  refers  to  the  case  where 
different  material  phases  exist  on  each  side  of  a  propagating  discontinuity. 
Moreover,  a  shock  wave  whose  Lagrangian  velocity  s  is  zero  is  referred  to  as  a 
contact  discontinuity  and  a  phase  boundary  with  s  =  0  is  called  a  stationary 
phase  boundary. 

Some  authors  define  shock  waves  as  traveling  discontinuities  that  satisfy 
the  so-called  subsonic-supersonic  condition  regardless  of  the  phases  on  each 
side  of  the  discontinuity.  The  subsonic-supersonic  condition  asserts  that  s  is 
bounded  below  by  the  isentropic  sound  speed  of  its  front  state  and  bounded 
above  by  the  isentropic  sound  speed  of  its  back  state  (see  Dunn  &  Fos- 
DICK  [11],  Truskinovsky  [27]).  For  a  material  with  a  nonconvex  Helmholtz 
energy,  it  is  possible  to  have  moving  discontinuities  where  the  phases  in  the 
front  and  the  back  are  different  but  the  speed  of  propagation  satisfies  the 
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subsonic-supersonic  condition.  Here,  though,  those  discontinuities  will  be 
referred  to  as  supersonic  phase  boundaries.  Truskinovsky  [27]  refers  to 
subsonic  phase  boundaries  as  kinks.  It  will  be  shown  that  only  subsonic 
phase  boundaries  have  room  for  a  kinetic  relation. 

2.3  Driving  traction  and  kinetic  relation. 

Using  (15)-(17)  in  (8)  provides  an  alternative  expression  for  the  energy  equa¬ 
tion,  i.e.. 


T]t  =  0  .  (24) 

It  follows  that  the  dissipation  inequality  (6)  is  trivially  satisfied  at  regular 
points.  The  rate  of  entropy  production  for  a  segment  [xi,X2]  of  the  bar 
which  contains  a  propagating  discontinuity  at  x  =  s{t)  can  be  expressed  as 
T{t)  =  Tbit)  +  Tsit),  where 

pX2 

^b{t)  =  /  fyqtdx  ,  Ts{t)  =  -pMs  .  (25) 

Jxi 

Here,  Tt  represents  the  bulk  entropy  production  and  Fg  corresponds  to  the 
entropy  production  due  to  the  moving  discontinuity.  Equations  (24)  and 
(25)  1  imply  that  F^  =  0,  hence  the  entropy  production  for  a  thermoelastic 
material  under  an  adiabatic  process  occurs  solely  because  of  the  presence  of 
a  moving  discontinuity.  Based  on  this  entropy  production,  define  the  driving 
traction  as 


f  =  -plvW)  ■ 


(26) 
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For  a  discussion  on  the  notion  of  driving  traction,  see  Abeyaratne 
Knowles  [1,  4]  and  Truskinovsky  [27].  For  a  thermoelastic  material, 
(14)-(17)  and  (26)  provide  the  following  equivalent  expression  for  the  driving 
traction  when  s  y^O 


/  =  P  {  M  -  (^7)  W  -  1^1 }  •  (27) 

In  terms  of  the  driving  traction,  the  dissipation  inequality  (12)  can  be  ex¬ 
pressed  as 


fs>0.  (28) 

Therefore,  the  jump  conditions  are  given  by  (9)-(12)  or,  equivalently,  by  (9), 
(13),  (14)  and  (28).  For  given  conditions  on  one  side  of  the  discontinuity,  say 
(7''',P'''),  and  for  a  given  speed  of  propagation  s,  the  set  of  points  (7“,^“) 
that  satisfy  the  jump  condition  (13)  is  called  the  Rayleigh  set.  Similarly,  for 
given  conditions  (7‘*‘,P'''),  the  set  of  points  (7“,^“)  that  satisfy  the  jump 
condition  (14)  is  called  the  Hugoniot  set  (see  Dunn  &  Fosdick  [11]).  The 
admissible  states  (7‘'',  6'^)  and  (7“,  0“)  are  those  that  belong  to  the  Rayleigh 
and  Hugoniot  sets  generated  at  either  (7''‘,^''')  or  (7“,^“)  and  also  satisfy 
the  entropy  jump  condition  (28).  In  problems  of  the  Riemann  type  for  bars 
made  of  thermoelastic  materials  that  can  undergo  phase  transformations,  the 
jump  conditions  are  not  sufficient  to  fully  determine  the  solution.  Additional 
information  is  provided  by  the  kinetic  relation.  Following  Onsager’s  approach 
to  irreversible  thermodynamics,  the  quantity  f/{0)  is  identified  as  an  affinity 
and  s  as  the  corresponding  flvx  (see  Abeyaratne  &  Knowles  [1,  4]).  The 


12 


kinetic  relation  is  part  of  the  constitutive  theory  and  relates  the  affinity  to 
the  flux  and  to  a  state  parameter,  in  this  case  the  temperature,  i.e., 

(29) 

where  F  is  a  materially- determined  function. 

2.4  A  specific  thermoelastic  material. 

In  order  to  obtain  results  in  specific  problems,  consider  a  thermoelastic  ma¬ 
terial  introduced  by  Abeyaratne  &  Knowles  [3].  A  phase  diagram  of 
this  material  in  the  7,0-plane  is  shown  in  Figure  1.  For  temperatures  below 
a  critical  temperature  6c,  the  material  can  exist  in  either  a  low  strain  phase 
Pi  or  a  high  strain  phase  P3.  These  phases  are  metastable  and  are  separated 
by  an  unstable  phase  P^.  Above  the  critical  temperature  the  material  can 
only  exist  in  a  stable  phase  P.  Throughout  this  analysis,  only  transforma¬ 
tions  from  or  to  the  low  and  high  strain  phases  are  considered.  A  thorough 
description  of  the  thermomechanical  characteristics  of  this  material  can  be 
found  in  [3].  The  boundaries  between  the  different  phases  are  given  by 

7m(^)  =  7c  +  M{e  -  Oc)  ,  7m(^)  =  7c  +  rn{d  -  do)  ,  (30) 
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where  7c  >  0,  0c  >  0)  and  m  are  constants.  The  expression  for  the 
Helmholtz  potential  is  given  in  each  phase  by 


Ay  -  —'y{0  -0t)-cO log  ( in  Pi, 
2p  p  \uT ) 


A 

2p 


.2  7r(7  -  IM? 

Tm 


V»(7, 0)  =  I 


——7(0  —  6t)  —  c01og 


OtJ 


in  P2, 


(31) 


A(^  -  7t)^  -  —(7  -  7r)(^  -  ^t) 

2p  p 

-c0  log  +  ^(0  -  6t)  in  P3, 


where  6t  is  the  transformation  temperature,  jt  is  the  transformation  strain, 
Xt  is  the  latent  heat  at  0  =  0x  and  7^  =  7m(^)i  7m  =  7m  (0)-  The  re¬ 
maining  parameters  were  defined  in  Section  2.2  and  are  assumed  constant. 
It  is  convenient  to  introduce  a  set  of  nondimensional  parameters  defined  as 
follows; 


T  = 


C0 


j.  7  -  V 

S  =  — ,  V  = - , 

7t  a7r 


_ 

~  ,2^2  > 


v  =  -,  M  = 
a 


a"7r 


G  =  Gjt, 


m  = 


(32) 


The  variable  6  might  be  viewed  as  a  “normalized”  strain.  The  parameter 
G  will  also  be  called  the  (modified)  Griineisen  coefficient.  Nondimensional 
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internal  energy  e,  stress  r  and  entropy  s  are  defined  as 


e  _  a  _  77 


(33) 


The  nondimensional  isentropic  sound  speed  is  ae  =  Ue/a,  hence,  from  (23), 
it  can  be  expressed  as  =  1  +  Q?T.  Other  nondimensional  parameters  will 
be  introduced  as  required.  From  (16)  and  (31),  the  stress  response  function 
is  given  by 


5-G{T-  Tt) 


f{S,T)  =  { 


S  —  Sm 


-  G(T  -  Tt) 


[(5-1)-G(T-Tt) 


low  strain  phase, 
unstable  phase, 

high  strain  phase. 


(34) 


where  =  ^m(T),  5m  =  5m  (T).  Observe  that,  for  a  fixed  temperature,  the 
stress-strain  relation  is  linear  in  each  phase.  Hence,  this  material  will  be 
referred  to  as  the  trilinear  material.  The  entropy  response  function,  from 
(17)  and  (31),  is  given  by 


r  /  T 

G5  +  f  1  -H  log  — 


low  strain  phase. 


s(5,r)  = 


G5  -1-  ^1  -I-  log 
1 

~2(m-  M) 


TtJ 

(5  -  Sc)^ 
[{T-TcY 


unstable  phase, 


(35) 


G(5  -  1)  +  ^1  -H  log  ^  high  strain  phase. 
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where  the  explicit  forms  of  5^  and  (5m  given  by  (30)  were  used.  The  jump 
conditions  for  the  special  thermoelastic  material  when  the  high  strain  phase 
is  on  the  right  of  the  phase  boundary  (and  the  low  strain  phase  on  its  left) 
are,  from  (13),  (14),  (28),  (31),  (34)  and  (35),  given  by 


(l-v")(i]-l-G(71=0,' 
|g{T)M  +  IT]  -  ({(S)  -  0  -  GTt  -  ir}  V  =  0  ,  , 

|g([(S1  -  l)  +  log^  -  ^1''  ^  0  ■ 


(36) 


Similar  expressions  can  be  obtained  when  the  low  strain  phase  is  on  the  right 
of  the  phase  boundary.  Using  (36)  i  in  (36)2  to  eliminate  [(5J  gives 


|■l(v)W  +  r2(v,T-)[Tl  +r3(v.r,T-)  =0  .  (37) 


where 


'■i(v) 

r2(v,T) 


G^v 


2(1  -  v2)  ’ 


r3(v,(5,T)  = 


V 

1  -V2 


(l-V^) 


-  +  GTt  +  It 


> 

(38) 


j 


The  dependence  of  r2  on  T  is  implicit  in  the  isentropic  sound  speed. 
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3  Augmented  theory 

3.1  Motivation. 


Although  the  classical  thermoelastic  theory  of  phase  transitions  has  room 
for — and,  in  fact,  requires — a  kinetic  relation,  such  a  piece  of  constitutive 
information  is  introduced  in  the  theory,  in  a  sense,  as  a  supplement.  Some 
other  theories  have  built-in  features  from  which  the  equivalent  of  a  kinetic 
relation  can  be  derived.  In  fact.  Fried  &  Gurtin  [14]  argue  that  a  kinetic 
relation  is  not  part  of  the  constitutive  theory  but  follows  as  a  consequence 
of  new  balance  principles.  Here,  nevertheless,  a  viewpoint  consistent  with 
Abeyaratne  &  Knowles  [1]  is  taken.  An  augmented  theory  that  includes 
viscosity,  strain  gradient  and  heat  conduction  is  introduced.  In  the  next 
section,  from  an  admissibility  criterion  for  a  special  class  of  solutions,  namely 
traveling  waves,  a  special  kinetic  relation  in  the  thermomechanical  case  is 
derived.  This  procedure  follows  well  established  techniques  used  in  fluid 
and  solid  mechanics  (see  Abeyaratne  &  Knowles  [2],  Hattori  [17], 
Rosakis  [22],  Slemrod  [24,  25],  Truskinovsky  [27]) 

By  using  the  classical  form  of  the  balance  principles  and  enforcing  the 
second  law  of  thermodynamics  for  all  admissible  processes,  Gurtin  [16], 
generalizing  a  result  established  by  Coleman  &  Noll  [10],  proved  that  the 
Helmholtz  potential  cannot  depend  on  higher  order  derivatives  of  the  strain. 
In  order  to  obtain  a  thermodynamically  consistent  theory  that  includes  strain 
gradient  and  viscosity  effects,  Dunn  &  Serrin  [12]  proposed  a  modification 
of  the  energy  balance  to  take  into  account  long  range  interactions  by  intro¬ 
ducing  the  concept  of  interstitial  working.  An  alternative  approach  is  used  by 
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Truskinovsky  [27].  In  this  section,  Dunn  &  Serrin’s  model  is  specialized 
for  the  one-dimensional  case. 


3.2  Basic  equations. 

In  the  regularized  theory,  one  introduces  from  the  onset  a  new  fundamental 
quantity  called  the  interstitial  working  p  by  DUNN  &  Serrin.  As  in  the 
classical  theory  of  Section  2,  the  momentum  equation  is  given  by  (1),  but  by 
introducing  the  interstitial  working  p  and  considering  heat  conduction,  the 
energy  equation  is 

+  P  +  =  p(^^  +  ^v^^dx  .  (39) 

The  Clausius-Duhem  inequality  (dissipation  inequality)  has  the  same  form 
as  in  the  classical  theory  with  heat  conduction,  i.e.. 

At  points  where  the  field  quantities  have  enough  differentiability,  the  lo¬ 
cal  form  of  the  balance  of  linear  momentum,  the  energy  equation  and  the 
Clausius-Duhem  inequality  are 


-  P'^t  5 


{av  +  q  +  p)x  =  p 


t 


(41) 

(42) 

(43) 
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and  the  compatibility  equation  is 


Vx  =  7t 


(44) 


Using  (41)  in  (42)  provides 


(TVx  +  gx+Px  =  P^t  ■  (45) 

Furthermore,  expanding  the  dissipation  inequality  (43)  and  using  the  expres¬ 
sion  of  Qx  provided  by  the  energy  equation  (45)  yields,  in  terms  of  the  free 
energy  potential  ip, 

P  {A  +  OtV)  -  (TVx-px-^  <0  .  (46) 

In  order  to  introduce  higher  gradient  effects  into  the  theory,  it  is  assumed 
that  Ip,  p,  T),  a,  q  are  given  by  functions  ip,  p,  fj,  a,  q  of  7,  7x,7xx)^i  and 

7t.  All  functionals  depend  on  the  same  set  of  variables  since  there  is  no  o 
priori  reason  to  discard  any  of  them.  Nevertheless,  some  restrictions  arise  in 
connection  with  the  second  law  of  thermodynamics. 

3.3  Thermodynamic  restrictions. 

The  previous  constitutive  assumption  is  compatible  with  thermodynamics  if 
every  admissible  thermodynamic  process  satisfies  the  dissipation  inequality 
(46).  This  requirement  restricts  the  functional  form  of  ip,  p,  f},  b  and  q. 
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Equations  (46)  and  (44)  provide 


’0'y'yt  d"  ('00  +  d“  d“  ‘^'txxl/xxt  d"  '00i^it  d"  '07t'7ttj 

q6 

iP'l'ix  d*  Pd^x  d"  P'ix'^xx  d"  P'yxxHxxx  d"  POx^xx  d"  P-ytl/xt)  ^  ^  0  • 


This  inequality  has  to  be  satisfied  for  all  deformations  and  all  temperatures 
at  any  point  {x,t).  It  is  possible  to  find  a  deformation  and  a  temperature 
field  for  which,  at  a  given  point  {x,t),  the  quantities  Oulxulxxu^xulu^lxxx 
and  6 XX  are  independent  of  the  values  of  7,^,7i,7xx,^®  and  7t  at  the  same 
point.  Therefore,  in  order  to  satisfy  the  dissipation  inequality,  the  following 
terms  must  vanish: 


00  d-  ^  =  0  ,  /907,  -  07,  =  0  ,  (47) 

07xx  =  0  .  0e.  =  0  ,  07,  =  0  ,  (48) 

07xx  =  0  ,  00.  =  0  .  (49) 

The  formal  role  of  the  interstitial  working  and  the  connection  between  the 
viscosity  and  the  strain  gradient  can  be  observed  from  (47)2-  Without  the 
interstitial  working  term  or  without  a  dependence  of  the  interstitial  working 
on  7t,  then  the  Helmholtz  potential  would  not  be  a  function  of  7i.  In  view 
of  (47)- (49),  the  dissipation  inequality  reduces  to 

—  \  X  f  \ 

P0^7t  -  a-7t  -  (077x  +  000X  +  07*7ix)  -  <  0  .  (50) 

Equation  (48)  implies  that  the  free  energy  cannot  depend  on  these  variables, 
hence  0  =  "0(7,  ^,7i)-  Similarly,  (49)  implies  that  p  —  0(7.^)7x,7t)-  Notice 
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that  (47)  1  corresponds  to  the  same  relation  between  entropy  and  free  energy 
as  in  the  classical  theory.  Since  •ip  does  not  depend  on  74,  then  (47)2  can  be 
integrated,  i.e., 


P(7,  7x,  7t)  =  (7,  7x)7t  +  P°(7,  7x)  ,  (51) 

where  p°  is  the  equilibrium  part  of  the  interstitial  work.  Using  this  expression 
for  the  interstitial  working  in  the  dissipation  inequality  (50)  gives 

(p  \i-i  -  ('^•r.)^]  -  7l  -  (p»  +  I)  ~  +P?.7.i)  <  0  . 

(52) 

The  third  term  in  this  inequality  is  independent  of  Ox  and  74;  therefore,  it 
should  satisfy  p®(7,^,7x)7x  +P°j7,^,7i)7ix  >  0,  for  all  7, 6»,  71,7^1.  Since 
')xx  can  be  specified  independently  of  7,  71  then,  necessarily,  =  0,  but 

then  p°  is  independent  of  7x  and  the  same  inequality  implies  p®  =  0.  There¬ 
fore,  p°  can  only  be  a  function  of  6  and  the  dissipation  inequality  (52)  becomes 

{p  [^7  -  7t  -  (pb  +  fj0x<0.  (53) 

The  equilibrium  values  of  a  and  q  can  be  obtained  by  replacing  9x  and  74 
by  hOx  and  /174  respectively  and  taking  the  limit  h  ->  0  in  the  dissipation 
inequality  (53),  i.e., 

{p  [^7  -  -  O-eq}  It  -  (Pe  +  -  ^ 
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where  Ceq  =  d-(7, 6, 7*,  0, 0)  and  q^q  =  9(7, 0, 7x,  0, 0).  Since  the  first  term  in 
(54)  is  linear  in  jt  and  the  second  is  linear  in  6x,  then  the  equilibrium  stress 
is  given  by 


^eq  P 


'07  {^Ix  ^  > 


and  the  equilibrium  heat  is 


(65) 


(56) 


For  elastic  materials  where  and  q  are  independent  of  7t,  the  first  term  of 
the  dissipation  inequality  (53)  is  linear  in  7*,  hence  the  elastic  part  of  the 
stress  is  completely  determined  by  the  free  energy  as 


tTei  =  p  (^7x)^  • 


(57) 


Notice  that  the  elastic  stress  is  equal  to  the  equilibrium  stress  (when  6x  =  0). 
Define  the  viscous  stress  as  cXy  =  a  —  cr^i  and  the  “non-equilibrium”  heat  flux 
as  ^ne  =  Q~Qeq-  fh©  general  case,  the  elastic  stress  can  be  obtained  from  the 
Helmholtz  free  energy,  but  for  the  viscous  stress  an  additional  constitutive 
assumption  is  required.  Returning  to  the  general  case,  using  the  expression 
(57)  for  the  elastic  stress  and  the  definitions  of  the  viscous  stress  and  the 
non-equilibrium  heat  flux,  the  dissipation  inequality  (53)  can  be  expressed 
as 


c^v7t  + 


Qne^x 

~e~ 


>  0 . 
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Moreover,  the  balance  of  linear  momentum  (41)  now  takes  the  form 


’ii  -  ]  +  (<''')•  -m  =  o 


(58) 


The  energy  equation  (45),  with  the  use  of  (15),  (47)  i  and  (51),  yields 


+  P  (^7x)  ^  7t  +  +  ?!  =  P^77t  +  P^Pt  ■ 

Prom  the  decomposition  of  the  stress  into  its  elastic  and  viscous  parts  and  the 
heat  flux  into  its  equilibrium  and  non-equilibrium  parts,  the  energy  equation 
now  becomes 


^v7t  (pfl)i  ^  (9ne)x  ~  P^Vt  ■  (59) 

This  is  the  classical  energy  equation  with  two  additional  terms:  the  work 
done  by  the  viscous  stress  (crv7t)  and  the  interstitial  working  term  (p^)^.  0- 
Moreover,  only  the  non-equilibrium  part  of  the  heat  flux  is  present. 

3.4  Specific  constitutive  assumptions. 

Introduce  the  following  constitutive  assumptions:  Assume  that  the  Helmholtz 
potential  is  of  the  form 

i’(7,  7x)  =  '0(7,  ^)  +  \^7l  ,  (60) 

where  A  is  the  strain  gradient  coefficient.  This  coefficient  is  formally  analo¬ 
gous  to  the  concept  of  capillarity  in  fluid  mechanics.  Moreover,  assume  that 
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the  viscous  stress  is  given  by 


(Tv  =  pu'Yt  ,  (61) 

where  u  is  the  viscosity,  and,  finally,  assume  that  the  equilibrium  heat  flux 
is  zero  and  the  non-equilibrium  heat  flux  is  only  a  linear  function  of  6x  (i.e., 
Fourier’s  law^) 

Q{^x)  -  9ne(^x)  =  kOx  ,  (62) 

where  k  is  the  thermal  conductivity.  For  simplicity,  the  strain  gradient  co¬ 
efficient,  viscosity  and  thermal  conductivity  are  assumed  constant  for  all 
phases.  Notice  that  the  viscous  stress  is  assumed  as  a  function  of  only  and 
the  non-equilibrium  heat  flux  as  a  function  of  dx  only.  These  assumptions 
are  consistent  with  the  thermodynamical  restrictions.  The  internal  energy 
of  the  regularized  theory  is  e  =  e(7,  d,  7^)  =  ^(7, 9, 7^)  —  (7, 9,  ^x)  which, 

by  (60),  can  be  expressed  as 

e  =  c(7, 9, 7*)  =  e(7, 9)  +  ^A7^  ,  (63) 

where  e(7, 9)  =  ^(7, 9)—9'ipe{'y, 9).  Moreover,  since  the  equilibrium  heat  flux 
is  zero,  then,  from  (56),  it  follows  that  the  equilibrium  part  of  the  interstitial 
work,  can  at  most  be  a  constant.  Since  the  energy  equation  and  the 
dissipation  inequality  only  involve  derivatives  of  then,  without  loss  of 

^Recall  that  in  the  linear  theory  (linearized  about  a  reference  temperature),  the  heat 
flux  can  only  depend  linearly  on  the  temperature  gradient  (see,  e.g.,  CARLSON  [7]). 
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generality  for  a  material  that  follows  Fourier’s  law,  one  can  take  =  0. 
With  assumptions  (60)-(62)  enforced,  the  stress  becomes 

d(7,  Tii,  It)  =  P  0)  -  A71X  +  nt]  ,  (64) 

the  balance  of  linear  momentum  (58)  takes  the  form 

{i’-y)  -  >^Jxxx  +  t'7xf  -Vt  =  0  ,  (65) 

and  the  energy  equation  (59)  reduces  to 

nt+^dxx  +  e{tP9)^  =  0.  (66) 

Recall  that  the  strain  and  particle  velocity  are  related  through  (44).  The 
dissipation  inequality  (53)  can  be  expressed  as 

pnl  +  ^  0  ,  (67) 

which  is  satisfied  if  1/  >  0  and  k  >  0.  By  setting  u  =  X  =  0  in  equa¬ 
tions  (65)  and  (66),  one  recovers  the  corresponding  equations  for  a  classical 
thermoelastic  material  with  heat  conduction. 
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4  Structured  wave:  Viscosity,  strain  gradient 


and  heat  conduction 

4.1  Generalities. 

The  traveling  wave  problem  is  analyzed  here  within  the  context  of  the  aug¬ 
mented  theory.  There  are  several  possibilities  for  connecting  the  augmented 
theory  of  Section  3  to  the  sharp  interface  theory  of  Section  2.  The  three 
structuring  mechanisms  are  related  to  viscous,  strain  gradient  and  heat  con¬ 
duction  eflFects.  The  case  when  the  viscosity  and  strain  gradient  coefficient 
are  assumed  small  (in  a  sense  to  be  specified)  is  anal3^ed  and  it  is  shown  that 
a  kinetic  relation  can  be  inherited  when  A,  i/,  fc  — >  0  and  that  the  augmented 
theory  reduces  to  the  sharp  interface  theory.  As  a  motivation  for  the  use  of 
viscous  and  strain  gradient  effects,  it  was  shown  in  [28]  that  when  the  only 
structuring  mechanism  is  heat  conductivity,  only  supersonic  traveling  waves 
can  be  achieved  and  no  kinetic  relation  can  be  derived.  It  is  also  important 
to  remark  that  the  presence  of  heat  conduction  induces  some  structure  to  the 
phase  boundary  but,  more  significantly,  it  permits  specification  of  different 
temperatures  at  -l-oo  and  —  oo.  This  in  turn  allows  an  appropriate  connection 
with  the  adiabatic  limit  case. 

After  introducing  some  general  features  of  a  structured  traveling  wave 
in  Section  4.2,  the  basic  equations  are  developed  in  Section  4.3.  Classical 
techniques  of  singular  perturbation  theory^  are  then  used  to  identify  the  re¬ 
gion  where  the  viscosity  and  strain  gradient  play  an  important  role  (i.e., 

^See,  e.g.,  Kevorkian  &  Cole  [18],  Whitham  [29]. 
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the  “inner”  solution)  and  the  region  where  heat  conduction  is  the  relevant 
structuring  mechanism  (i.e.,  the  “outer”  solution).  A  uniformly  valid  ap¬ 
proximation  is  then  obtained  from  the  asymptotic  matching  of  the  inner  and 
outer  solutions  and  a  kinetic  relation  of  the  form  (29)  is  derived  for  a  subsonic 
traveling  phase  boundary. 

4.2  Traveling  wave:  Preliminaries. 

A  steady,  structured  traveling  wave  corresponds  to  the  case  when  all  fields 
satisfy  (41)-(44)  and  depend  only  on  the  variable  ^  =  x  -  st,  where  s,  the 
wave  speed,  is  a  constant.  A  moving  discontinuity  can  be  viewed,  within  the 
context  of  the  adiabatic,  sharp  interface  theory,  as  the  limit  of  a  structured 
wave  within  the  framework  of  the  augmented  theory.  Hence,  in  an  infinitely 
long  bar,  the  conditions  at  ±oo  for  a  structured  wave  in  the  augmented  theory 
should  correspond  to  the  front  and  back  states  of  a  moving  discontinuity  in 
the  theory  without  augmentation.  The  precise  statement  of  the  traveling 
wave  problem  will  be  given  in  Section  4.4.  In  this  section,  some  general 
characteristics  of  the  traveling  wave  are  recorded  for  future  use.  Notice  that 
although  in  the  classical  theory  of  Section  2  the  function  -ip  corresponds  to 
the  Helmholtz  free  energy,  in  the  augmented  theory,  the  free  energy  is  given 
by  •0.  Assume  there  exists  a  structured,  steady  traveling  wave  for  ^  €  E. 
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From  (41),  (42)  and  (44)  it  follows  that 


a'  +  psv'  =  0  ,  (68) 

u'  +  57'  =  0  ,  (69) 

{ov)'  +  q'  +p'  +  ps  ^e  +  =  0  ,  (70) 

where  (•)'  denotes  differentiation  with  respect  to  Integration  of  these 
equations  yields 


cr  +  psv  =  Hi  ,  (71) 

v  +  s'y  =  H2  ,  (72) 

av  +  q+p  +  ps(^  +  =  Hs  ,  (73) 

where  the  Hi  axe  integration  constants.  Let  the  conditions  at  ^  >  ±00  be 

given  by 


7(-oo)  =  7  , 

7(+oo)  =  7+  , 

(74) 

?;(— 00)  =  v~  , 

u(+oo)  =  v'^  , 

(75) 

1 

II 

1 

${+00)  =  6'^  , 

(76) 
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and  all  derivatives  of  these  quantities  vanish  when  ^  — >  ±oo.  Hence,  the 
integration  constants  must  satisfy 


Hz  =  cr+t;+ 


Hi  =  <T^  -{■ 

H2  =  v'^  s'y'^ 

-\-ps  + 


=  a  +  psv  , 

(77) 

=  v~  +  57“  , 

(78) 

=  a~v~  +  ps 

> 

(79) 

where  =  <7(7^,  0"^,  0,0)  and  -  9^'4>0{'y^,6^,O).  With 

the  foregoing  constitutive  assumptions  (60)-(62),  the  stress  cr"^,  given  by  (64), 
is  <7^  =  pip^{'y^,6^)  and  the  internal  energy  e*  becomes  =  ■^(7^,0^)  — 
Therefore,  from  (77)-(79),  the  values  given  in  (74)-(76)  must 
satisfy  the  adiabatic  jump  conditions  (9)-(ll).If  a  traveling  wave  solution 
exists  then,  from  (59),  it  necessarily  satisfies 

7/2  QH 

-psT]'  =  • 


Integrate  from  ooto^— ^•ooto  get 

where  [r/]  =  r]+  -  t]~  =  -•0e(7'^,^'*')  +  ‘^b{T ,0-)  .  Integrating  by  parts  the 
second  term  in  the  integral  and  using  the  conditions  at  ±00  provides 


To  make  the  connection  with  the  classical  adiabatic  theory,  consider  the  driv¬ 
ing  traction  /  given  by  (26).  Thus,  the  values  in  (74)-(76)  must  also  satisfy 
fs>0,  which  corresponds  to  the  entropy  jump  inequality  (28).  The  previ¬ 
ous  analysis  shows  that  the  field  values  at  ^  ±oo  for  the  traveling  wave 
problem  must  satisfy  all  the  corresponding  jump  conditions  of  the  classical 
adiabatic  problem.  Moreover,  the  condition  for  the  existence  of  a  traveling 
wave  solution  for  the  augmented  theory  should  provide  an  additional  relation 
for  the  phase  boundary  velocity  s  in  terms  of  the  driving  traction  /. 

4.3  Basic  equations. 

Evaluate  the  constants  Hi  in  (71)-(73)  using  the  conditions  at  an  arbitrary 
point  ^  =  ^0-  These  equations  can  be  written  as  follows: 


a  -  (To  +  ps{v  —  Vo)  =  0  ,  (80) 

{v  -  Vo)  -I-  s(7  -  7o)  =  0  ,  (81) 

{crv  -  o-qUo)  +q-qo  +  p-Po  +  ps(€-€o  +  ^{v^  -  vo^)^  =  0  , 

(82) 

where  (Tq  =  ^(7o,^o,7o>  ~^7o)>  and  70  =  7(^0))  etc.  The  velocity  u  can 
be  eliminated  from  (80)  and  (82)  by  use  of  (81).  The  momentum  equation 
(80)  becomes 


(T  -  ps'^'y  =  cro-  psSo  , 


(83) 
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and  the  energy  equation  (82)  can  be  expressed  as 


q-qo  +  p-po  =  -s 


p{€  -  Co)  -  (  ^  t'  —  )  (7  -  7o) 


(84) 


Alternatively,  using  the  expression  for  a  given  by  (83),  the  previous  equation 
can  be  written  as 


q-qo+p-po  +  ps 


(e  -  eo)  -  Y  (7  -  lof  -  y (7  -  7o) 


=  0  . 


(85) 


Equation  (85)  is  therefore  a  combination  of  the  compatibility,  momentum 
and  energy  equations.  Moreover,  it  is  clear  that  by  choosing  the  conditions 
behind  (resp.  ahead)  a  right-traveling  wave  in  (77)-(79),  i.e.,  the  ”  side 
(resp.  the  “-t-”  side),  (83)  and  (85)  are  also  valid  with  the  subscript  0  replaced 
by  ”  (resp.  by  “-1-”)  and  with  q~  =  p~  =  0  (resp.  g"*"  =  p"*"  =  0).  Using 
the  constitutive  assumptions  (60)-(62),  the  heat  flux  q  and  the  interstitial 
work  p  are,  for  a  traveling  wave,  given  by 

q  =  kO'  ,  p  =  —psX'y'^  .  (86) 

Prom  (64)  and  (63),  the  stress  and  the  internal  energy  are  given  by 


a  =  p  [^.y(7, 9)  -  A7"  -  vsi] 
e  =  e(7,0)-t-^A7'^  • 


(87) 

(88) 
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Using  (86)-(88)  in  the  governing  equations  (83)  and  (85)  gives 


Kl"  -  7o)  +  -  7o)  -  [^7(7,  ^)  -  V'7(7o,  ^0)]  +  s^(7  -  7o)  =  0  , 

(89) 


-9'q)  +  s  |e(7, e)  -  e(7o, ^o)  -  ^7'^  -  7^^) 

-\s^  (7  -  laf  -  [^7(70,  ^0)  -  A70  -  1^57^  (7  -  7o) I  =  0  .  (90) 

Introduce  a  new  function  w  =  7'.  Equations  (89)  and  (90)  can  be  expressed 
as  a  system  of  first  order  ordinary  differential  equations: 

f 

i  =  w, 

w'  -Wo  =  -J  {us{w  -  Wo)  -  [^^(7,  e)  -  •0^(70, 9o)]  +  s^(7  -  7o)}  , 

< 

•  r  A  -2 

9'  -9'o  =  |e(7, 9)  -  e(7o,  9o)  -  -(w^  -  “  y  (7  -  lof 

-  [^7(70, 9o)  -  Xw'o  -  i/swo]  (7  -  7o) }  • 

The  following  stability  analysis,  with  some  modifications,  follows  from  Ngan 
&  Truskinovsky  [20].  Let  (7*,tt;*,0»)  be  an  equilibrium  point  of  this  dy¬ 
namical  system.  Therefore,  from  (91)  and  assuming  s  ^  0,  these  points 
satisfy  the  jump  conditions  of  the  classical  adiabatic  theory  for  a  given  state 
{'y'^,9'^)]  see  relations  (13)  and  (14).  The  states  ai  ^  ±00  given  by  (74)- 

(76)  are  equilibrium  points.  To  analyze  the  nature  of  the  equilibria,  consider 
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the  linearized  system  around  an  equilibrium  point,  i.e.,  W  =  AW,  where 
W  —  {'y,w,6)  and 


A  = 


us 

T 


A 


k 


(92) 


where  tp*g  =  'ipyei'y*,^*),  Voe  =  Notice  that  to 

obtain  the  term  Az\,  the  relation  e-y  =  ip-y  —  Oip^Q  was  used.  The  term  .^33 
was  obtained  using  eg  =  —O’ipeg.  Let  l3  =  Pi,  i  =  1,2,3,  he  the  eigenvalues  of 
A.  The  characteristic  polynomial  of  >1  is  —  hP"^  +  I2P  —  /3  =  0,  where  li, 
z  =  1, 2, 3,  are  the  principal  invariants  of  >1.  A  direct  calculation  gives,  using 
the  material  parameters  introduced  in  Section  2.2,  gives 


Define 


H  =  h  + 


hh 

3 


Lz  —  i\-\-  4^2  , 


(96) 


then,  when  tz  <  0,  all  three  eigenvalues  are  real  and  different;  when  4.3  =  0,  all 
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eigenvalues  are  real  and  two  of  them  are  equal  and  when  ta  >  0  one  solution  is 
real  and  the  two  other  are  complex.  The  eigenvalues  A  satisfy  /^i+A+Z^a  =  h 
and  A/^2/^3  =  h,  hence  if  the  three  eigenvalues  are  real,  and  since  /i  <  0 
when  s  >  0,  then  the  three  eigenvalues  cannot  be  all  positive.  Therefore,  if 
h  >  0,  the  equilibrium  saddle  point  has  a  two-dimensional  stable  manifold 
and  a  one-dimensional  unstable  manifold.  If  s  <  0,  then  the  equilibrium 
point  has  a  two-dimensional  unstable  manifold  and  a  one-dimensional  stable 
manifold.  If  /a  <  0  then  the  situation  is  reversed. 

4.4  Traveling  phase  boundary. 

Since  the  fields  under  consideration  are  piecewise  differentiable,  the  traveling 
phase  boundary  problem  has  to  be  divided  into  three  sub-problems.  Suppose 
that  the  conditions  at  ^  >  +00  correspond  to  the  low  strain  phase  and  the 
conditions  at  ^  >  —00  correspond  to  the  high  strain  phase.  Presumably  there 
is  an  interval,  say  (0,6),  where  the  points  are  in  the  unstable  phase.  Prom 
the  above  considerations  and  in  view  of  the  smoothness  requirements  of  the 
traveling  wave,  the  solution  must  satisfy  conditions  (74),  (76)  and  the  strain, 
strain  gradient  and  temperature  must  be  continuous  at  ^  =  0  and  (  =  b. 
In  particular,  7(0=*")  =  7m(^o)  and  7(6=*=)  =  7m(^6),  where  the  functions  ‘jrn 
and  7a/  are  given  by  (30)  and  the  temperatmes  6q  =  ^(0^)  and  9b  =  ^(6^) 
are  unknown.  Additionally,  for  the  class  of  traveling  waves  considered  here, 
the  phase  segregation  conditions  must  be  enforced,  i.e.,  the  solution  must  be 
such  that  the  material  is  in  its  high  strain  phase  in  (— oo,0),  unstable  phase 
in  (0,6)  and  low  strain  phase  in  (6, -l-oo).  Moreover,  the  temperature  has  to 
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be  below  the  critical  value  9c.  These  conditions  can  be  expressed  as 


-1  <  7(0  <  7m(9(0  for  -  00  <  ^  <  0  , 
7m(^(0)  <  7(0  <  7m(9(0)  for  0  <  ^  <  6  , 

7(0  >  7mW0)  for  6  <  ^  <  +00  , 

^(0  ^  for  —  00  <  ^  <  +00  . 


(97) 


Structured  phase  boundary  traveling  wave  problem:  Given  one  end 
state,  say  at^  —00,  (or,  alternatively,  {'y'^,6'^)  at  ^  ^  +00)  find 

all  admissible  states  {'y'^,6'^)  at  ^  +00  in  the  high  strain  phase  such  that 

{‘y~,9~)  and  (7"'’,^''’)  can  be  connected  via  a  traveling  wave.  A  connection 
between  the  end  states  is  achieved  if  there  exist  functions  7(^)  and  9(^), 
^  G  R,  that  satisfy  the  phase  segregation  conditions  (97)  and  the  following 
three  problems: 


Problem  1:  Low  strain  phase.  Suppose  that  for  ^  G  (—00, 0)  the  material 
is  in  its  low  strain  phase.  Using  the  conditions  at  —00  in  system  (91)  and  the 
corresponding  expressions  of  ^(7, 0),  c{'y,6),  e~  and  a~  given  by  equations 
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(15)-(17)  and  (34),  the  governing  equations  are 


y  =  w, 

w'  =  — Y  {I'sw  —  (a^  —  s^)(7  —  7“)  +  aa^{9  —  0")}  , 
A 

O'  =  +  aa^r  (7  -  7") 

+0(6  -  $-)  - 

subject  to  the  following  boundary  conditions: 

7(-oo)  =  7"  ,  ^(-00)  =  ,  7'(-oo)  =  0  , 

7(0")  =  7o  =  7c  +  M{6q  -  6c)  ,  [^0  =  ^(0")]. 

The  temperature  60  is  not  known  a  priori. 


Problem  2:  Unstable  phase.  Suppose  that  for  ^  G  (0,  b)  the  material  is  in 
its  unstable  phase.  The  governing  equations  are 

( i  =  'u}. 


^ vsw  —  {a?  —  s^)(7  —  7  )  +  aa^{9  —  6) 

2  7r(7-7M(^))] 

'  7m(^)  -  7m{0)  }  ’ 

(101) 


9'  =  ~  )  +  c{9-9  )  - 

])■ 


{29  -  9c)  {'y  -  ycf 


2{M  -  m)  [  {9-  9c)^ 


—  M^9c  —  2M(7  —  'yc) 


subject  to  the  following  initial  conditions: 


7(0+)  =  70  =  7c  +  M(0o  -  9c)  ,  0(0+)  =  0o  ,  (102) 


y(o+)  =  7'(o-) . 


(103) 
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Problem  3;  High  strain  phase.  Suppose  that  for  ^  G  (6,  +oo)  the  material 
is  in  its  high  strain  phase.  The  governing  equations  axe 

y  =  w, 

w'  =  — i  {usw  —  {a?  —  s^)(7  —  7'*')  +  (xa?{6  —  , 

A 

ff  =  ~  ~  ('•'  “  - 1*) 

+c{e-e*)-^w^Y 

subject  to  the  following  boundary  conditions: 

7(6+)  =  75  =  7^7  +  m{eb  -  dc)  ,  [^6  = 

7(+oo)  =  ,  ^(+00)  =  6'^  ,  7'(+oo)  =  0  , 

y(6+)  =  7'(6-)  =  7^  . 

The  temperature  db,  the  strain  gradient  y  and  the  width  of  the  unstable 
phase  b  are  unknown  a  priori. 

4.5  Nondimensional  parameters. 

It  is  convenient  to  introduce  a  length  parameter  based  on  the  viscosity  and 
strain  gradient  coefficient  and  to  use  it  to  obtain  nondimensional  equations. 
Recall  that  the  viscosity  is  given  in  terms  of  L^T~^,  where  L  is  length  and 


(105) 

(106) 
(107) 
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T  is  time,  and  the  strain  gradient  coefficient  in  terms  of  L^T  Let 

pea  V  k 

Therefore,  the  system  (98)  can  be  expressed  as 
5'  =  w, 

<  w'  =  —uwMSN  +  [(1  —  v^)((5  —  5~)  —  G(T  —  T~)]  , 

(108) 

r"  =  -V  {i  (1  -  V^)  («  -  S-f  +  GT-{S  -  S-)  +  {T-T-)  -  , 

where  (•)'  =  d/dz.  The  boundary  condition  at  z  =  0  is 

S{0)  =  6c  +  M{T{0)-Tc)  .  (109) 

The  nondimensional  expression  for  the  unstable  phase  will  be  given  in  sec¬ 
tion  4.6.  The  system  (104)  can  be  written  as 

5'  =  w, 

^  w'  =  —uw\/w  +  -aP'  [(1  —  v^)(5  —  (J'*')  —  G(T  —  T"*")]  , 

(110) 

r  =  -V  |i  (1  -  v^)  (i  -  +  GT*(i  -  r)  +  (r-T*)  -  I , 
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and  the  boundary  condition  (105)  at  z  =  b  is 


5{h)  =  5c  +  m{T{h)-Tc)  .  (Ill) 

The  boundary  condition  (107)  will  be  analyzed  in  Section  4.7.  The  dimen¬ 
sionless  driving  traction  is  defined  as 

hence,  if  the  conditions  at  -foo  correspond  to  the  high  strain  phase,  then, 
from  (35), 

f=  -  G(P]-l)  +  logP-|;  (T).  (112) 

Prom  (108)2,  the  temperature  can  be  written  in  terms  of  the  strain  and  its 
derivatives  as 

r(z)  -  T-  =  i  [{1  -  v^)({(z)  -  r)  -  la,vJ'(z)  -  ^r(z)]  , 

(113) 

hence,  equation  (108)3  can  be  expressed  as  a  single  equation  for  the  strain, 
i.e., 

r  -  li<5"  +  12^'  -  h{5  -S-)  +  ^  [(<5')'  -  -  v"*)  (<5  -  S-f]  =  0  . 

(114) 

A  similar  equation  can  be  obtained  for  the  high  strain  phase. 
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4.6  Perturbation  analysis. 

To  obtain  an  approximation  of  the  solution  in  the  metastable  phases,  consider 
a  singular  perturbation  of  the  equations.  Let 


1 


and  consider  the  case  where  A  <C  1  ,  -C  1  but  the  ratio  a;  remains  fixed. 
Hence,  assume  that  e  is  a  small  parameter.  For  definiteness,  the  analysis  is 
carried  out  for  the  low  strain  phase  (it  is  formally  the  same  for  the  high  strain 
phase).  As  a  general  scheme  of  notation  in  this  section,  an  over  bar  represents 
a  quantity  related  to  the  inner  expansion  whereas  a  quantity  without  an 
over  bar  usually  represents  a  value  related  to  the  outer  expansion.  A  prime 
denotes  differentiation  with  respect  to  the  outer  coordinate  and  a  dot  denotes 
differentiation  with  respect  to  the  inner  coordinate. 

Outer  solution 

Assume  that  the  asymptotic  expansion  for  6  corresponds  to  a  power  series  in 
s  for  fixed  z  ^  0,  i.e., 

5(z)  —  (5“  ~  0o(z)  +  +  •  •  •  as  e  — >•  0  . 

This  corresponds  to  the  outer  solution.  Using  this  representation  for  5  in 
equation  (114),  collecting  terms  in  powers  of  e  and  setting  each  of  the  coef¬ 
ficients  equal  to  zero  gives,  for  the  first  term, 

(p'o  =  r20o  +  ^o4'o  ) 
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where  r2  is  defined  by  (38)  and 


roCv)  =  .  (115) 

There  are  two  functions  that  satisfy  the  equation  for  (one  of  them  being  the 
trivial  solution).  For  the  nontrivial  solution  the  first  term  of  the  expansion 
is  given  by 

•^0(2)  =  -  |-coexp[-r2z]  -  l|  ,  (116) 

I'D  L  •'0  ) 

where  co  is  a  constant.  In  order  to  enforce  the  condition  that  0o(z)  ->  0 
as  z  —00,  it  is  required  that  r2  =  r2(v,  T")  >  0  in  the  low  strain  phase, 
but  no  further  restriction  is  placed  upon  Cq.  For  the  high  strain  phase  the 
restriction  is  that  r2  =  r2(v,  T"^)  <  0.  Observe  that  Vn  6  N,  (t>^\z)  — >  0  as 
z  —00.  The  same  occurs  for  z  — >■  +00.  Formally,  the  strain  gradient  is 
given  by  5'{z)  ~  ^o(z)  +  . . . ,  and  the  temperature  is  given,  from  (113),  by 

T(z)-T-~l^A,(z)  +  ...  ;  (117) 

therefore,  the  boundary  condition  T(z)  —  T“  -»  0  as  z  —00  is  satisfied. 

Inner  solution. 

For  the  inner  solution  (close  to  z  =  0),  consider  a  new  coordinate  defined  as 

z 

x  =  -  =wz  . 
e 
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Assume  the  following  asymptotic  expansion  for  fibced  x: 


5(x)  —  5  ~  <^o(x)  +  £</>i(x)  +  . . .  . 

Proceeding  as  in  the  previous  section,  the  equation  to  determine  the  first 
term  of  the  expansion  is,  from  (114), 


00  +  va;0o  -  (1  -  v^)0o  =  0  • 

The  solution  0o  is  given  by  the  sum  of  two  exponentials  plus  a  constant.  The 

i 

exponents  are  /3x,  where 

P  =  i(-va;  ±  Vv2a;2  +  4(1  -  v^))  .  (118) 

The  exponents  can  be  either  real  or  complex.  For  the  subsonic  range  (sub¬ 
sonic  with  respect  to  the  isothermal  sound  speed),  the  exponents  are  real: 
y3>0,  /3<0for0<v<l  and  /3  >  0  for  -1  <  v  <  0.  In  the  low  strain 
phase  (where  x  <  0),  since  the  inner  and  outer  solution  must  match  in  some 
region  (as  explained  below),  the  exponential(s)  with  negative  exponent(s) — 
which  would  provide  an  unbounded  term  in  the  limit — should  be  discarded. 
For  the  high  strain  phase,  the  solution  with  positive  exponent(s)  should  be 
abandoned.  If  both  exponents  have  the  appropriate  sign,  then,  in  general, 

0o(x)  =  coe^"  +  cie^"  +  c2  .  (119) 

Notice  that,  from  (96)  in  nondimensional  form,  the  eigenvalues  of  the  system 
linearized  about  the  equilibrium  point  for  the  metastable  phases  are  real  and 
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different  if 


-^(1  -  [vV  +  4(1  -  v^)]  +  0(e)  <  0  , 

i.e.,  if  v^a;^  +  4(1  —  v^)  >  0  for  small  e.  Therefore,  for  small  e,  the  equilib¬ 
rium  points  of  the  systems  (108)  (low  strain  phase)  and  (110)  (high  strain 
phase)  have  a  two-dimensional  stable  manifold  and  a  one-dimensional  un¬ 
stable  manifold  for  0  <  v  <  1  and  a  one-dimensional  stable  manifold  and  a 
two-dimensional  unstable  manifold  for  — 1  <  v  <  0. 

Unstable  phase. 

For  the  unstable  phase,  the  system  of  equations  (101)  in  nondimensional  form 
and  expressed  in  inner  coordinates  becomes 

f  j  =  w. 


W  =  — (JVW  -t- 


(i-v2)(5-r)-G(T-r-)  + 


5-Sc-M{T-Tc) 

{M-m){T-Tc) 


r  =  -^  (j  (1  -  v")  (« -  i-?  +  Gr-(i  -  S-)  +  {T-  T-)  - 


(120) 


2(M-m)L  {T-Tcy 


1} 


Assuming  that  the  strain,  strain  gradient  and  the  temperature  are  bounded 
throughout  the  unstable  phase,  it  is  possible  to  conclude  from  (120)3  that 


r(x)  -  To  =  0{e)  , 
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where  To  is  the  (unknown)  value  of  the  temperature  at  x 
assuming  that 

5(x)  -  5"  ~  0o(x)  +  £^i(x)  +  . . .  , 
equations  (120)i,2  give,  for  the  first  term, 

^o(x)  +  ^i0o(x)  +  A:2  =  0  ) 


where 


h 


(l-v^)  + 


1 

(M-m)(ro-Tc) 


A:2  =  G(To-T-)- 


S- -5c-U{T^-Tc) 

(M-m)(ro-rc) 


Hence, 


<^o(x) 


=  c5‘e^“='  +  ^ 


where 


/3„  =  ^(-vo;  ±  y/v^uj^-4ki)  . 

Notice  that  the  exponents  can  be  either  real  or  complex, 
quantities  will  be  useful  in  the  foregoing  analysis: 

P(v)  =  13  —  0  =  +  4(1  —  v^)  , 

Q(v,ro)=^„-^„  =  Vv2a;2-4/ci  .  ^ 


0.  Therefore, 


(121) 

(122) 

(123) 
The  following 

(124) 
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4.7  Subsonic  case. 


To  obtain  a  uniformly  valid  approximation  of  the  solution  in  the  metastable 
phases,  the  inner  and  outer  solutions  must  match.  For  the  subsonic  case 
(subsonic  with  respect  to  the  isothermal  sound  speed)  the  regions  where  the 
different  approximations  are  valid  are  shown  in  Figure  2. 

Low  strain  phase. 

Suppose  that  0  <  v  <  1.  From  (88)2  it  follows  that  r2(v,  T*)  <  0,  hence 
the  nontrivial  outer  solution  for  the  low  strain  phase  has  to  be  abandoned  in 
favor  of  the  trivial  solution.  Furthermore,  since  /3  >  0  and  /3  <  0,  the  inner 
solution  for  the  low  strain  phase  involves,  up  to  0(s),  only  the  exponential 
with  exponent  /3x.  From  (119),  the  inner  solution  is  given  by 

5(x)  -S~  =  ^o(x)  +  0(e)  ,  (125) 

where  4>o{x)  =  +  C2.  Before  enforcing  the  boundary  condition  at  x  =  0, 

it  is  convenient  first  to  analyze  the  restrictions  imposed  by  the  asymptotic 
matching  of  the  inner  and  outer  solutions.  The  inner  expansion,  rewritten 
in  outer  variables  (i.e.,  z)  and  expanded  for  small  e  and  fixed  z  is  given  by 

C2  + _  The  one-term  outer  expansion  of  the  one-term  inner  expansion 

rewritten  in  inner  coordinates  is,  clearly,  C2.  For  this  expansion  to  match  the 
(trivial)  outer  expansion,  it  is  required  that 

C2  =  0  .  (126) 
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The  temperature  can  be  expressed,  from  (113)  in  inner  coordinates,  as 

r(x)  -  T"  =  i  (1  -  v^)((5(x)  -  S~)  -  u>vS(x)  -  5(x)j  .  (127) 

The  strain  and  the  temperature  at  x  =  0  are  related  through  the  bound¬ 
ary  condition  (109).  Prom  (125)-(127),  the  first  term  of  the  temperature’s 
expansion  at  x  =  0  is 

^  (1  -  v2)0o(O)  -  a;vJo(0)  “  ?o(0)]  =  ^(1  “  =  0  , 

and  the  first  term  of  the  strain’s  expansion  at  x  —  0  is  hence,  from  (109) 
and  matching  coefficients  of  the  same  order  gives,  for  the  first  term, 

^  =  —S  -|-  (5c  +  M(T  —  Tc)  .  (128) 

Therefore,  the  strain,  strain  gradient  and  the  temperature  at  z  =  0  are  related 
to  the  conditions  at  z  — >•  —  oo  by 

^(0)  —  S  =  Co  +  0(e)  , 

m  =  0co  +  O(e)  ,  [  (129) 

T(0)  -  T- =  0(e)  ,  ^ 

± 

where,  for  given  (5~,T~)  and  v,  cq  is  given  by  (128)  and  (3  by  (118).  Observe 
that  the  strain  gradient  is  given  here  with  respect  to  the  inner  coordinate. 
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In  view  of  the  above,  the  composite  expansion  for  the  low  strain  phase  is 


5{z)  —5  =  +  0(e)  for  z  <  0  . 


Therefore,  for  given  {6  ,T  )  and  v,  the  composite  expansion  for  the  low 
strain  phase  is  completely  determined  up  to  0(e). 

High  strain  phase. 

4* 

For  the  inner  solution  close  to  z  =  b,  since  /?  >  0,  the  exponential  term  with 

-f 

exponent  fix  is  discarded.  Hence,  the  inner  solution  is 


(5(x)  -  J"'’  =  (f>o{x)  +  0(e)  , 


where  ^o{x)  =  +  C2  and  b  =  b/e.  The  boundary  condition  (111)  at 

X  =  b  relates  the  strain  and  the  temperature.  Prom  (127),  expressed  for  the 
high  strain  phase,  it  follows  that  the  lowest  order  term  of  the  temperature’s 
expansion  at  x  =  b  is 


1 

G 


(1  -  v^)0o(b)  -  a)v0„{b)  -  (j||(b)  =  1(1  -  V^)C2  . 


(130) 


The  first  term  of  the  strain’s  expansion  at  x  =  b  is  Ci  +C2,  hence,  using  (130) 
the  lowest  order  term  in  (111)  gives 


Therefore,  the  strain,  strain  gradient  and  the  temperature  at  x  =  b  are  related 
to  the  conditions  at  z  — )■  +oo  by 


5(b)  —  5'*'  —  Cl  +  C2  +  0(s)  , 
i(b)  =  ^c,  +  0(£), 
r(B)-T+  =  i(l-v2)c2  +  0(£). 


(132) 


In  this  case,  the  outer  solution  is  nontrivial  and  the  first  term  is  given  by 
(116).  For  the  high  strain  phase  it  is  convenient  to  work  with  z  —  b  instead 
of  z.  The  translation  is  absorbed  by  the  integration  constant  cq.  The  outer 
expansion,  rewritten  in  inner  variables  and  expanded  for  small  e  (with  x 
fixed)  gives 

— - - +  ...  , 

r2Co  -  ro 

where  r2  =  r2(v,  T"*").  The  inner  expansion  rewritten  in  outer  coordinates, 
expanded  for  small  e  (with  z  fixed)  and  rewritten  in  inner  coordinates  gives 
C2  H-  . . . ;  hence,  for  the  inner  and  outer  expansion  to  match  asymptotically, 
it  is  necessary  that 


r2 

C2  =  -  . 

r2Co  -  ro 

The  one-term  composite  expansion  is  obtained  by  adding  the  one-term  inner 
and  outer  expansions  and  subtracting  the  one-term  outer  expansion  of  the 
one-term  inner  expansion,  i.e.. 
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The  temperature  in  the  high  strain  phase  is,  from  (117),  (127)  and  the  above 
expression  for  the  strain,  given  by 


r(z)-T+  =  -(l-v2) 

VD 


r2 


r2Coe-''=*(^-‘’)  -  roj 


+  0{e) 


Notice  that  the  inner  expansion  of  the  temperature  is,  up  to  0(e),  a  constant 
given  by  (1/G)(1  -  v^)c2. 


Unstable  phase. 

Since  To  =  T(0~)  =  T“  +0(e),  then,  from  (122)  and  enforcing  the  continuity 
of  the  strain,  strain  gradient  and  temperature  at  x  =  0,  it  follows  that 


q  +  ^--jf-=Co  , 
ki 

. 


Let 


and  set  =  5m(T  )  and  5m  =  5m{T  );  using  the  identities 


5m  ~  ^ 


h 

kk 


(134) 
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it  follows  that 


c?  =  i 


Q 


^iPu-P) 


Co  =  - 


Q 


+ 

-  k 

0 


(^M  ~  ^u)  , 


(135) 


where  Q  =  Q(v,r-)  is  given  by  (124),  P  by  (118),  by  (134)  and  (121) 

db  _ 

and  P^  by  (123).  The  phase  segregation  condition  at  x  =  b  specifies  that  the 
material  is  on  the  interface  between  the  unstable  and  the  high  strain  phase. 
Since  the  temperature  throughout  the  unstable  phase  is  essentially  constant 
(up  to  0{e)  )  then,  from  (122)  and  up  to  the  leading  term,  this  condition 
corresponds  to 


^  =  s^-Su. 


(136) 


Now,  enforcing  continuity  of  temperature  at  x  =  b,  gives,  from  (132)3, 


Using  the  above  expression  in  (131),  the  constant  ci  becomes 
Cl  =  Sc-S^  +  m{T-  -  Tc)  -  -T+)  . 

Let 

K  =  (137) 

1  —  v^ 
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This  quantity  can  be  thought  of  as  the  strain  (in  the  high  strain  phase)  cor¬ 
responding  to  the  front  state  of  an  isothermal  phase  boundary  for  a 

given  back  state  (5~,T~).  The  state  {5^,T~)  can  also  be  viewed  as  the  back 
state  of  an  adiabatic  shock  wave  for  a  given  front  state  It  is  im¬ 

portant  to  notice  that  the  strain  gradient  at  (5;f ,  T~)  does  not  vanish,  hence 
this  state  is  not  an  equilibrium  point  of  the  system  (110).  Using  (36)i  and 
(137),  the  constants  ci^  can  be  expressed  as  ci  =  5m  —  and  C2  =  5^  —  5'^. 
The  boundary  condition  (111)  was  used  to  determine  the  asymptotic  expan¬ 
sion  in  the  high  strain  phase  and  (136)  was  obtained  by  enforcing  the  phase 
segregation  condition  for  the  asymptotic  expansion  in  the  unstable  phase. 
Therefore,  the  continuity  of  the  strain  at  x  =  b  is  automatically  satisfied 
from  the  continuity  of  the  temperature.  The  only  remaining  condition  that 
has  not  been  enforced  is  the  continuity  of  the  strain  gradient  at  x  =  b  (i.e., 
(107)).  Prom  the  inner  expansion  for  the  strain  in  the  high  strain  phase,  it 
follows  that 

5(b+)  =  Pci  +  0{e)  =  ;9(5m  -  St)  +  0{£)  . 

On  the  other  hand,  from  the  asymptotic  expansion  in  the  unstable  phase, 
one  has 

i(b-)  =  +  0(e)  ; 

hence,  up  to  the  leading  term,  the  closing  condition  for  the  construction  of 
a  traveling  wave  is 

-  K)  ■  (138) 
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It  is  possible  to  prove  that 


-  Su 


u 


In  view  of  the  above  relation  and  using  (124)  and  (135),  equations  (136)  and 
(138)  can  be  expressed,  after  some  simplifications,  as 


and 


g-VO-b  ^ 


(140) 


Since  <  1,  then  P  >  0.  Thus,  if  Q  >  0,  equation  (139)  has  no  solution 
since  the  left-hand  side  is  always  greater  than  the  right-hand  side  for  any 
b  >  0.  Therefore,  a  necessary  condition  to  obtain  a  solution  is  Q  =  i\Q\-  Of 
all  the  possible  roots  of  (139),  only  the  smallest  satisfies  the  phase  segregation 
conditions,  hence,  choosing  the  argument  of  the  complex  number  P4-i|(5|  to 
be  in  the  interval  (0, 7r/2),  then,  from  (139), 


b  =  i^[27r-4arg(P-+z|Q|)]  .  (141) 

The  width  of  the  unstable  region  is  therefore  determined  as  a  function  of  the 
velocity  and  the  temperature  T~ . 

The  phase  space  for  the  system  of  equations  (108),  (110)  and  (101)  is 
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(5,  5',  T).  Figure  3  shows  the  projection  into  the  6-T  plane  of  the  trajec¬ 
tory  of  a  traveling  wave  such  that  T'^  >  T~  (in  this  case  v  =  0.05).  The 
Hugoniot  curve  T  =  T^{S),  as  defined  in  Section  2.3,  is  also  shown.  The 
corresponding  wave  profiles  are  shown  in  Figures  4  (inner  coordinates)  and  5 
(outer  coordinates).  These  graphs  were  obtained  using  the  following  values: 
G  =  0.8,  M  =  1.4,  m  =  1.2,  Sc  =  3.1,  Tx  =  1,  Tc  =  3,  lx  —  0.5,  u)  =  1  and 
£  =  \/w  =  10"®.  Similarly,  Figures  6,  7  and  8  show  a  traveling  wave  such 
that  T'*'  <  T"  (in  this  case  v  =  0.7).  Observe  that,  in  this  case,  the  first 
portion  of  the  wave,  from  J"*"  to  5^,  is  expansive. 

Kinetic  relation. 

As  mentioned  in  Section  4.2,  the  values  (5*,T*,v)  are  formally  related  via 
the  jump  conditions  (36)i,2.  If  equation  (141)  is  used  to  compute  b,  then 
(140)  can  be  viewed  as  an  additional  restriction  on  the  values  of  (J*,r*,v). 
Hence,  for  a  given  pair  of  values  (e.g.,  back  state  {S~,T~)  or  front  state 
or  other  combinations),  the  system  (36) i, 2,  (141)  provides  the  nec¬ 
essary  restrictions  to  determine  the  other  three  values  (e.g.,  ((J"*", T"*",  v),  etc.). 
Notice  that  this  relation  only  holds  for  0  <  v  <  1  and,  from  the  entropy  in¬ 
equality,  it  is  also  required  that  the  driving  traction  should  be  positive.  The 
condition  (140)  can  also  be  interpreted  in  terms  of  the  driving  traction.  To 
this  end,  the  strain  S^  can  be  written,  from  (140),  as 


1  -  (/5//5)  exp[-va;b/2] 
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Therefore,  from  (121),  (134)  and  (137),  the  strains  6  and  can  be  ex¬ 
pressed  as 


6- 


1  /  ~  \ 
(1_V2) 


+  ) 


and 


5t  = 


+ 


(1  -  V2)  V^M  -«5n 

Hence,  formally,  since  b  is  a  function  of  v  and  T~  and  in,  and  5m  are  functions 
of  T~,  then  5^  =  5u(v,  T~)  and  6~  =  5"(v,  T~).  Using  this  expression  for  6~ 
in  (37),  the  temperature  T"*"  can  be  obtained  as  a  function  of  v  and  T~.  The 
strain  5''‘  in  the  high  strain  phase  can  be  determined  from  (36)  i.  Finally, 
from  (112),  it  is  possible  to  obtain  a  kinetic  relation  of  the  form  (29),  i.e.. 


f  =  f(v,T-). 


It  is  also  possible  to  use  T"*"  as  a  reference  temperature  by  solving  for  T~ 
and  5^  for  given  T'^  in  which  case  f  =  f(v,  T"^).  Several  kinetic  curves,  cor¬ 
responding  to  different  values  of  oj  and  a  common  value  of  T'^  =  1.5  are 
shown  in  Figure  9.  The  last  point  of  each  curve  (for  maximum  v)  corre¬ 
sponds  to  the  limit  5+  =  5ni  as  required  by  the  phase  segregation  condition 
(97)3.  All  curves  have  a  common  value  at  v  =  0  (the  parameter  u  is  always 
multiplied  by  v  in  the  equations).  Observe  that  for  small  enough  values  of 
u),  the  corresponding  kinetic  curve  is  non  monotonic.  Relations  of  this  kind 
have  recently  been  considered  by  Rosakis  &  Knowles  [23]  in  a  purely 
mechanical  context.  It  is  interesting  to  point  out  that  setting  v  =  0  in  (108), 
(110)  and  (120)  results  in  T(z)  being  constant,  which  is  different  from  the 
limit  V  — >•  0  in  the  solution  obtained  by  perturbation  analysis.  Moreover,  the 
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energy  jump  condition  (86)2  and  the  entropy  inequality  (86)3  are  trivially 
satisfied  when  v  =  0.  Figure  10  shows  a  set  of  kinetic  relations  for  different 
values  of  T'*'  and  a  common  value  of  a;  =  1.  In  this  case,  for  a  given  velocity 
V,  the  driving  traction  decreases  for  increasing  temperature  T'^ — the  same 
behavior  is  observed  when  T~  is  used  as  a  reference  temperature. 

The  negative  velocity  case  is  different  than  the  positive  one  in  the  sense 
that  /?  and  /3  have  the  same  sign.  Even  though  it  is  possible  to  obtain  a 
uniformly  valid  expansion  for  the  low  strain  phase,  there  seems  to  be  no 
simple  way  to  construct  such  an  expansion  for  the  high  strain  phase. 

An  interesting  point  that  arises  from  this  analysis  is  that,  although  it 
is  possible  to  formally  cast  the  admissibility  criterion  in  terms  of  a  kinetic 
relation  of  the  type  proposed  in  [5],  there  is  no  simple  dependence  on  the  tem¬ 
perature.  The  temperatme  enters  the  kinetic  relation  in  a  somehow  arbitrary 
way  in  the  sense  that  either  the  temperature  in  front  of  a  phase  boundary  or 
the  temperature  in  the  back  can  be  used  as  a  reference  temperature.  This 
fact  was  also  observed  by  Ngan  h  Truskinovsky  [20]  in  their  numerical 
simulations  and  suggests  that  the  functional  form  of  the  kinetic  relation  in 
the  thermoelastic  case  should  be  reconsidered. 


References 


[1]  R.  Abeyaratne  and  J.  K.  Knowles.  On  the  driving  traction  acting  on  a 
surface  of  strain  discontinuity  in  a  continuum.  Journal  of  the  Mechanics 
and  Physics  of  Solids,  38(3):345-360,  1990. 

[2]  R.  Abeyaratne  and  J.  K.  Knowles.  Implications  of  viscosity  and  strain- 
gradient  effects  for  the  kinetics  of  propagating  phase  boundaries  in 
solids.  SIAM  Journal  on  Applied  Mathematics,  51(5):1205-1221,  1991. 

[3]  R.  Abeyaratne  and  J.  K.  Knowles.  A  continuum  model  of  a  thermoelastic 
solid  capable  of  undergoing  phase  transitions.  Journal  of  the  Mechanics 
and  Physics  of  Solids,  41(3):541-571,  1993. 

[4]  Rohan  Abeyaratne  and  J.  K.  Knowles.  Kinetic  relations  and  the  prop¬ 
agation  of  phase  boundaries  in  solids.  Archive  for  Rational  Mechanics 
and  Analysis,  114:119-154,  1991. 

[5]  Rohan  Abeyaratne  and  J.  K.  Knowles.  Dynamics  of  propagating  phase 
boundaries:  Adiabatic  theory  for  thermoelastic  solids.  Physica  D, 
79:269-288,  1994. 

[6]  H.  B.  Callen.  Thermodynamics  and  an  Introduction  to  Thermostatistics. 
John  Wiley,  2“*^  edition,  1985. 

[7]  D.  Carlson.  The  linear  theory  of  thermoelasticity.  In  S.  Fliigge,  editor, 
Handbuch  der  Physik,  VI  2/a.  Springer- Verlag,  1972. 

[8]  J.  W.  Christian.  The  Theory  of  Transformations  in  Metals  and  Alloys. 
Pergamon  Press,  2"*^  edition,  1981. 

[9]  R.  J.  Clifton.  On  failure  waves  in  glasses.  Applied  Mechanics  Reviews, 
48:540-546,  1993. 

[10]  B.  D.  Coleman  and  W.  Noll.  The  thermodynamics  of  elastic  materials 
with  heat  conduction  and  viscosity.  Archive  for  Rational  Mechanics  and 
Analysis,  13:167-178,  1963. 

[11]  J.  E.  Dunn  and  R.  Fosdick.  Steady,  structured  shock  waves.  Part  1: 
Thermoelastic  materials.  Archive  for  Rational  Mechanics  and  Analysis, 
104:295-365,  1988. 


57 


[12]  J.  E.  Dunn  and  J.  Serrin.  On  the  thermodynamics  of  interstitial  working. 
Archive  for  Rational  Mechanics  and  Analysis,  88:95-133,  1985. 

[13]  J.  L.  Ericksen.  Introduction  to  the  Thermodynamics  of  Solids.  Chapman 
and  Hall,  1991. 

[14]  E.  Pried  and  G.  Grach.  Rate  and  gradient  effects  in  a  theory  for  solid- 
solid  phase  transition.  In  K.  R.  Rajagopal,  editor.  Recent  advances  in 
elasticity,  viscoelasticity  and  inelasticity.  World-Scientific,  1995. 

[15]  E.  Pried  and  M.  Gurtin.  Dynamic  solid-solid  transitions  with  phase 
characterized  by  an  order  parameter.  Physica  D,  72:287-308,  1994. 

[16]  M.  Gurtin.  Thermodynamics  and  the  possibility  of  spatial  interaction  in 
elastic  materials.  Archive  for  Rational  Mechanics  and  Analysis,  19:339- 
352,  1965. 

[17]  H.  Hattori.  The  Riemann  problem  for  a  van  der  Waals  fluid  with  en¬ 
tropy  rate  admissibility  criterion.  Isothermal  case.  Archive  for  Rational 
Mechanics  and  Analysis,  92:247-263,  1986. 

[18]  J.  Kevorkian  and  J.  D.  Cole.  Multiple  Scale  and  Singular  Perturbation 
Methods.  Springer- Verlag,  1996. 

[19]  M.  Lusk.  Martensitic  phase  transitions  with  surface  effects.  Journal  of 
Elasticity,  34:191-227,  1994. 

[20]  S-C.  Ngan  and  L.  Truskinovsky.  Kinetic  relations  for  adiabatic  phase 
boundaries.  In  preparation. 

[21]  O.  Penrose  and  C.  Pife.  Thermodynamically  consistent  models  of  phase- 
field  tjqje  for  the  kinetics  of  phase  transition.  Physica  D,  43:44-62,  1990. 

[22]  P.  Rosakis.  An  equal-area  rule  for  dissipative  kinetics  of  propagat¬ 
ing  strain  discontinuities.  SIAM  Journal  on  Applied  Mathematics, 
55(1):100-123,  1995. 

[23]  P.  Rosakis  and  J.  K.  Knowles.  Unstable  kinetic  relations  and  the  dy¬ 
namics  of  solid-solid  phase  transitions.  In  preparation. 


58 


[24]  M.  Slemrod.  Dynamics  of  first  order  transitions.  In  M.  Gurtin,  editor, 
Phase  Transformations  and  Material  Instabilities  in  Solids.  Academic 
Press,  1984. 

[25]  M.  Slemrod.  A  limiting  “viscosity”  approach  to  the  Riemann  problem 
for  materials  exhibiting  change  of  phase.  Archive  for  Rational  Mechanics 
and  Analysis,  105:327-365,  1989. 

[26]  C.  Truesdell.  Rational  Thermodynamics.  Springer- Verlag,  2"^*  edition, 
1984. 

[27]  L.  Truskinovsky.  Kinks  versus  shocks.  In  R.  Fosdick  J.E.  Dunn  and 
M.  Slemrod,  editors.  Shock  Induced  Transitions  and  Phase  Structures  in 
General  Media.  Springer- Verlag,  1991. 

[28]  S.  Turteltaub.  Dynamics  of  Phase  Transformations  in  Thermoelastic 
Solids.  PhD  thesis,  California  Institute  of  Technology,  1997. 

[29]  G.  B.  Whitham.  Linear  and  Nonlinear  Waves.  John  Wiley,  1974. 


59 


Low  strain  phase  Unstable  phase  High  strain  phase 

Figure  1:  Material  phases  in  the  temperature-strain  plane. 
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Figure  2:  Outer  and  inner  solutions  for  metastable  phases. 
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Figure  3:  Projection  of  the  trajectory  of  a  traveling  wave  from  the  6, 5',  T- 
phase  space  into  the  S,  T-plane  for  T"*"  >  T  .  The  Hugoniot  curve  T  =  (6) 

is  based  at 
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Figure  6:  Projection  of  the  trajectory  of  a  traveling  wave  from  the  5, 5',  T- 
phase  space  into  the  S,  T-plane  for  T'*’  <T~.  The  Hugoniot  curve  T  =  {5) 

is  based  at 
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T,d 


Figure  7:  Strain  and  temperature  wave  profiles  in  inner  coordinates  for  v  = 
0.7. 
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T,b 


Figure  8:  Strain  and  temperature  wave  profiles  in  outer  coordinates  for  v  = 
0.7. 
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